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A quantum analogue of the Jarzynski equality is constructed. This equality connects an en- 
semble average of exponentiated work with the Helmholtz free-energy difference in a nonequi- 
librium switching process subject to a thermal heat bath. To confirm its validity in a practical 
situation, we also investigate an open quantum system that is a spin 1/2 system with a scanning 
magnetic field interacting with a thermal heat bath. As a result, we find that the quantum 
analogue functions well. 

KEYWORDS: quantum nonequilibrium dynamics, Jarzynski equality 



Understanding nonequilibrium quantum dynamics well be an important problem in the next 
century; new types of quantum devices, for instance, may function in their nonequilibrium states. 
In fact, it has been theoretically shown that quantum ratchet systems, which have the potential 
for being a new type of quantum device, function well in competitive situations between thermal 
relaxation and external driving, which are extremely nonequilibrium conditions J3i'i) A framework 
of theoretical methods for understanding such systems, however, is still under development. In 
this letter, we shall develop one equality which holds during a nonequilibrium process accompanied 
with thermal relaxation and external driving. This equality is a quantum analogue of the Jarzynski 
equality Considering such equality raises important issues regarding quantum nonequilibrium 
thermodynamics. 

Consider the situation in which a classical system interacts with an isothermal heat bath char- 
acterized by an inverse temperature (3. In this case, the system will become a thermal equilibrium 
state after a sufficiently long relaxation time. Next we assume that the system depends on a 
(macroscopic) parameter A. If we change the parameter A infinitely slowly from to 1, the system 
becomes another thermal equilibrium state corresponding to A = 1 from an equilibrium state with 
A = 0. When the operation is sufficiently slow (switching time ts is much larger), that is, when the 
system remains in quasistatic equilibrium during a switching process, the total work W{ts = oo) 
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performed on the system by the outside through the parameter A is identical to the Helmholtz 
free-energy difference AF between the initial state and the final one: 



W{ts = oo) = AF = Fo - Fi , 

where F\ is the Helmholtz free energy of the thermal equilibrium state with A. If we operate the 
switching much faster, that is, a finite ts, the work must be larger than the free-energy difference: 

W{ts) > AF . 

This inequality represents the least- work principle. Extra work is required for a faster switching 
process because of thermodynamic stability. This extra work is released into the heat bath after 
the operation, when the system interacts with the isothermal heat bath. 

In 1997, Jarzynski presented an equality which differs from the above inequality in the same 
situation.i'& His result is expressed by the following form: 

{exp{-(3W{t,)))path = exp(-/3AF) . (1) 

This is called the Jarzynski equality. The brackets (. . .)path mean taking the ensemble average along 
a single path during the switching process. An initial state is chosen from a thermal equilibrium 
ensemble with /3 and A = 0. The dynamics of the system during the process can be taken as 
an arbitrary one that can express the thermal relaxation. If Monte Carlo dynamics, for instance, 
is used, one path corresponds to a single sampling series. Jarzynski also showed that the above 
equality holds in several dynamics such as Monte Carlo dynamics, Nose-Hoover dynamics, Langevin 
dynamics, and even Hamiltonian dynamics without a thermal heat bath.& This equality is general 
since it involves the least-work principle mentioned above. 

In this letter, we consider an analogue of the Jarzynski equality in quantum situations. Let us 
define several notations and quantities. We denote the Hamiltonian of the system as 'Hx(t), and 
A(t) a switching parameter operated by the outside at time t. A partition function Zx and the 
Helmholtz free energy F\ are expressed as follows: 

=Trexp(-/37^;,(i)) , (2) 
Fxit) =-p-' In Zxit){P) . (3) 

The work operator W{t) during a time interval [0,t] is defined by the following: 







ds d\ 

If A is not a macroscopic parameter, this definition does not correspond to the intuitive view of 
work. Here, we only consider the case where A is a macroscopic parameter. Using a density matrix 
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of the system p{t), we get the expression of the total average work {W{t)) in the switching process 
as 

dX{s) dHx(s) 



In the present situation, the least-work principle also holds 

{W{t)) > AF 

for arbitrary switching processes. The equal sign is valid for quasistatic limits as well as the classical 
case. 

Let us consider a path average of exponentiated work corresponding to the left-hand side of the 
classical Jarzynski equality. In the present situation, taking the path average is difficult because 
we do not know how an ensemble reproduces thermal relaxation in the quantum dynamics. There- 
fore, we use the density matrix approach here. Then, the classical path average can be directly 
interpreted into the quantum situation: 



N-l 



(exp {-(5W{t)))path « exp (-(iW{t)) ^ Tr lim T \{ {v%^^^^e-^'^'^^-+^'> e'^^'^^^^A 

(6) 

where bt = tjlSi and t„ = nbt. The term e ^'^^'•^n+i) q+I^'^^m corresponds to exponentiated in- 
finitesimal work. In this expression, "P^j^) is an infinitesimal time-evolution superoperator for a tiny 
step 6t at time t defined by 

-pSt _ C{t)5t 

where C{t) is the Liouville superoperator governing the time evolution of the density matrix: 

Details of the dynamics are not important at this time. Its least property is determined in the 
following. T shows that the time ordering is taken as an increase in the left direction. 
This is actually the analogue of the classical path average. It can be expressed as follows: 

{eM-PWmpath = lim dx^Yll dx„Pf(,^^^)(x„+i|x„)e-^^-(-") , (7) 

where x„ denotes canonical variables describing the classical system and an infinitesimal work, 
bw(xn) is defined by 'HA(t„+i)(2;n) ~'^\(tn){^n)- T^x\t _^_-^)i^n+i\xn) is an infinitesimal time-evolution 
operator such as a transition probability for the Monte Carlo dynamics. From this form we can 
understand that eq. @ is a straight extension of classical situations. 

In this way, we obtain the quantum analogue of the path-averaged exponentiated work. In the 
next step, we must show that the above expression is identical to the exponentiated Helmholtz 
free-energy difference. For this purpose, we use the fundamental property of the Liouville super- 
operator: it is vanished by the one density matrix that is the thermal equilibrium distribution 



3 



oc exp (^—PHx^t)^, because the system always interacts with the thermal heat bath in the present 
case. Such distribution, of course, depends on time t. That is, exp (^—/3Tix(^t)) is a singular solution 
of the dynamics. In this sense, a part of the exponentiated infinitesimal work does not evolve: 

A(tn+i; 

Finally, we obtain 

TV Jim T n {P%^^,)e-^'^^(*r.+,)e.+P-H^i.tr.)^^ 



n=0 ^A(O) 



,,(0) Ijv-oo ^(*^) J 
—TV I lim Vfif se"''^^(*iv)e+''^^(*^-i) ...Pfl >e-^^M*2)e+^^M*i)e-^^M*i)l 



^A(O) 



Zx(o) 



— Tr le-'^^^wj 

^A(O) 

Thus, the following equality 



exp 



(-l3W{t 



Zxio) 



is proofed. This is the quantum analogue of the classical Jarzynski equality, the "quantum Jarzynski 
equality." 

It should be noted that this expression is valid in pure quantum dynamics as in classical dynamics. 
In pure dynamics, an operation of the supcroperator on an arbitrary operator A, +i)^' 
identical to e ^'^^^^n+i)^^/^'- j^g'^Htn+i)^^/^ _ Then, all the infinitesimal time-evolution superopcrators 
are effectively identity operators. Finally, we obtain the exponentiated frcc-cncrgy difference again. 

The quantum Jarzynski equality is formulated by an operator form. Thus, all the variables 
during the process are operators in contrast to the classical one. This property creates some 
difficulties in the decomposition of the time variable. If we change the definition of the infinitesimal 
exponentiated work slightly, for example, to e'^^'^^'^^"'>e ^'^^'•^ri+i) or e ^^{*n))^ causes 

errors at each time. Denoting the error of order 0{5t) as At„, that is, the exponentiated work is 
taken as e"^^^(*"+i^e+^^^(*") + At„ + 0{{dtf), we obtain the final error as 

( N~l / JV-1 * ^ \ ^ 



' -Tr 



^A(O) 



1=0 \ n=l+l 



(8) 

where Il^z+i* represents the omission of the term Il^Jjv- limit N ^ oo this error may 

remain, since N terms of the order 0{5t) = 0{1/N) are collected. The higher-order errors may 
vanish because their number is 0{N) at most. Therefore, the small mistaken choice of the exponen- 
tiated infinitesimal work causes a serious error in the final result. In a practical situation, however, 
Af^ may be of the order 0{5t x (ts)~^)- For sufficiently large tg, the final error is not serious. 
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To confirm the quantum Jarzynski equahty, we numericahy investigate a simple quantum system. 
We choose the spin 1/2 system interacting with an isothermal heat bath and time-dependent 
magnetic field. Here, we consider a process in which the magnetic field is linearly reversed in a 
finite time interval. The Hamiltonian T~Ls{t) and the magnetic field A(t) at time t are as follows: 

-Hsit) = -^a^ - \Aa^ , 
X{t) = Ao (^1 - 2^ 

where and are Pauli matrices. The parameters A and Aq are constants taken as A = 0.1 
and Aq = 1, respectively. Interaction between the system and the heat bath is described by 
70"^ Z]o(oJt+aa) (7 is a coupling constant) and the bath is expressed by a set of harmonic oscillators 
J2a (^aOa + 1/2^ , where a\^{aa) is a creation (annihilation) operator with a mode a. tg is a 
switching interval. The dynamics starts at time t = with a thermal equilibrium state and finishes 
at time t = ts- Using the projection techniqueil* with the assumption that the heat bath is always 
in a thermal equilibrium state with an inverse temperature /?, we obtain the equation of motion of 
the density matrix of the spin system p{t) as 



^^=-^[nsit),p{t)] + ^rip) , (9) 



where F is a relaxation term guaranteeing that the system will relax into an instantaneous thermal 
equilibrium state oc exp{—(37is{t)). If the Hamiltonian does not depend on time, this equation 
can express the thermal relaxation to a thermal equilibrium state. The explicit form of T is so 
complicated that is not presented here, but can be essentially described by (interaction between 
the spin and the bath) and the autocorrelation function of the thermal heat bath variables. For 
actual calculations we use the second-order perturbation for the coupling constant 7. 

Figure 1 shows the results of numerical calculations. Several time series of the expectation value of 
work defined by eq. (|5[) are presented. Each time series has a different switching time t^. The exact 
Helmholtz free-energy difference at the corresponding time and the result of the quantum Jarzynski 
equality are also shown. As switching time increases, the value of work is more converged into the 
exact Helmholtz free-energy difference, because the dynamics becomes closer to a quasistatic one. 
On the other hand, when tg is short, a nonadiabatic transition dominates the dynamics. Thus, the 
extra work is needed. The result of the quantum Jarzynski equality completely agrees with the 
exact Helmholtz free-energy difference. 

To summarize, we have constructed a quantum analogue of the Jarzynski equality. This equality 
connects a type of average of the exponentiated work operator during the switching process with 
the exact Helmholtz free-energy difference between an initial thermal equilibrium state and the 
final thermal equilibrium state, even though the actual final state of the process is not thermal 
equilibrium. We have also confirmed that the present equality works in practical calculations. In 



5 




Fig. 1. Time series of the expectation value of work, the exact Helmhohz free-energy difference and the result of 
the quantum Jarzynski equality. Horizontal axis is taken to be scaled time t/ts- 



the spin 1 /2 system with varying magnetic field interacting with the thermal heat bath, the result 
coincides with the exact difference of the Helmholtz free energy. 

The present quantum Jarzynski equality includes the classical quality. In the classical limit, the 
density operator becomes a diagonal matrix, that is, the distribution function of the corresponding 
classical system. In this limit, the infinitesimal time-evolution operator describes only the transition 
between diagonal elements such as one of a Fokker-Planck operator. Then, the average part of the 
quantum one is similar to eq. (0). 

It should be noted that the dynamics during the switching process can be chosen arbitrarily as 
one which has only the property that the time-dependent thermal equilibrium state oc eyi'p{— jSTixit)) 
is a singular solution of the dynamics. This property does not restrict the dynamics to Liouvil- 
lian dynamics with the projection technique such as eq. (P). Even though the singular solution 
exp(— /?7^;^(t)) is unstable, the dynamics having such a singular solution can produce the same 
results. In this situation, however, other results derived from the dynamics are not physical. 

It is known that the classical Jarzynski equality is related to the fluctuation theorem,Bll& which 
is the relation of a distribution function of the entropy production rate, and is also valid in dy- 
namics between nonequilibrium steady states .0 Based on these findings, we can put forth several 
questions. What is a quantum analogue of the classical fluctuation theorem and how is the quan- 
tum Jarzynski equality extended into nonequilibrium steady dynamics? These problems are now 
in progress. 
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